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Motivated by spin-wave continuum (SWC) observed in recent neutron scattering experiments 
in Herbertsmithite, we use Gutzwiller-projected wave functions to study dynamic spin structure 
factor ^(q, cj) of spin liquid states on the kagome lattice. Spin-1 excited states in spin liquids 
are represented by Gutzwiller-projected two-spinon excited wave functions. We investigate three 
different spin liquid candidates, spinon Fermi-surface spin liquid (FSL), Dirac spin liquid (DSL) and 
random-flux spin liquid (RSL). FSL and RSL have low energy peaks in ^(q,^;) at K points in the 
extended magnetic Brillouin zone, in contrast to experiments where low energy peaks are found at 
M points. There is no obviuos contradiction between DSL and neutron scattering measurements. 

Besides a fractionalized spin [i.e. spin-1/2), spinons in DSL carry a fractionalized crystal momentum 
which is potentially detectable in SWC in the neutron scattering measurements. 


Quantum spin liquid states has been catching more 
and more attention in condensed matter physics[l-3]. 
They are new states of matters that are beyond the de¬ 
scription of Landau’s symmetry breaking theory of con¬ 
ventional ordered phases [4]. Spin degrees of freedom 
in quantum spin liquids are not frozen at zero tem¬ 
peratures, but highly entangled with one another over 
long ranges. The symmetries in long-range entangled 
many-body systems can be fractionalized[5, 6]. Quantum 
spin liquids allow deconfined spinon excitations which 
carry a fractional spin (z.e. spin-1/2) and give rise to 
spin-wave continuum (SWC) through a pair of spinon 
particle-hole excitations[7-10]. In some spin liquid states, 
spinons carry “fractional crystal momenta”. As a re¬ 
sult, the momentum resolved density of states for spin- 
1 excitation continuum {i.e. two-spinon excitations) 
has a period smaller than the elementary Brillouin zone 
(BZ) which is potentially detectable in neutron scattering 
measurements)!, 6]. 

Herbertsmithite [ZnCu 3 (OH) 6 Cl 2 ], a layered spin-1/2 
kagome lattice antiferromagnet, is a promising com¬ 
pound for an experimental realization of spin liquid 
states [11-20]. Recent inelastic neutron scattering on sin¬ 
gle crystals of Herbertsmithite[19] has detected a dif¬ 
fuse low energy SWC over a large energy and mo¬ 
mentum regions. Nearest-neighbor kagome antiferro¬ 
magnetic Heisenberg model (KAFHM) has been sug¬ 
gested for spin-liquid physics observed in Herbert¬ 
smithite. Many different ground states have been 
proposed for KAFHM [21-35]. Recent density-matrix- 
renormalization-group calculations [36-38] support a Z 2 
gapped spin liquid ground state and indicate the kagome 
ground state is proximate to a critical state. The low- 
energy spin excitations are also studied for these pro¬ 
posed candidate ground states[26, 29, 39-43] and in the 
exact diagonalization[44]. 

In this letter, we will compute dynamic spin structure 
factor S'(q, w) for spin liquid states on the kagome lattice. 
The ground state for a spin liquid is described by the 


Gutzwiller-projected wave function (GPWF) by project¬ 
ing out double occupancy components in the mean field 
ground state [1, 45]. Similarly, GPWFs for spin-1 excited 
states are constructed by applying Gutzwiller projection 
onto spinon-antispinon excited wave functions[46, 47]. As 
well as equal-time spin factor 5'(q) in the ground-state 
GPWF, we use Monte Carlo method to calculate the pro¬ 
jected Hamiltonian system {H, O} where IHI and O are the 
Hamiltonian matrix and wave function overlap matrix, 
respectively, in a subspace consisting of spin-1 spin-wave 
excited states[46, 47]. The projected Hamiltonian system 
is diagonalized through the general eigen equation which 
gives eigenvalues as the spinon-antispinon excitation en¬ 
ergies and spectrum representation for S'(q, w)[46, 47]. 

The best variational GPWF for the ground state of 
KAFHM is the Dirac spin liquid (DSL) [24, 28, 35]. DSL 
has flux TT in the hexagons of kagome lattice in the mean 
field Hamiltonian. For comparison, we also study a zero- 
flux state which is a spinon Fermi-surface spin liquid 
(FSL). If the spin system in Herbertsmithite doesn’t 
reach a true ground state, a random-flux spin liquid 
(RSL) is also possible. We find that all three spin liq¬ 
uid states have a SWC spectrum in S'(q, w) with a low 
intensity in the elementary BZ and high intensity in 2nd 
BZ. The spectrum width of SWC is around ~ 3 J and the 
integrated intensity of S'(q, w) up to 0.2J corresponds to 
around 20% of the equal-time spin structure factor S (q). 
Unlike one-dimensional (ID) antiferromagnetic spin-I/2 
chain[48], the bottom boundary edge of SWC is weakly 
dispersive and the intensity at the edge SWC is not di¬ 
vergent. Above the low boundary edge, S'(q, w) is almost 
energy independent and weakly depends on the momen¬ 
tum over a wide range of momentum. These general fea¬ 
tures of SWC in spin liquids on the kagome lattice are 
consistent with experimental observations. FSL has a 
low energy gap in S'(q, w) at the M point and low energy 
intensity peaks at K points in the magnetic BZ (MBZ). 
RSL has a similar S'(q, w) to FSL, but the gap at M 
points is smeared out. 
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In comparison, the low energy intensity peaks of 
«S'(q, w) in the experiments[19] are located at M points in 
the MBZ and high intensity region connecting M points 
goes through M" points instead of K points. Therefore, 
FSL and RSL are not likely to be spin liquid states real¬ 
ized in Herbertsmithite. DSL has no obvious conflict in 
^(q, w) with neutron scattering measurements. Particu¬ 
larly, the momentum resolved density of states for spin-1 
excitations in DSL has two cones at low energies at M 
and M” points. Different from FSL, spinons in DSL carry 
a fractionalized crystal momentum. As a result, the mo¬ 
mentum resolved density of states for spin-1 excitations 
is periodic in one-quarter (shadow region in Fig. 1 (d)) of 
the elementary BZ. The boundary edge below SWC for 
DSL resembles the mean field continuum edge and Dirac 
cones around M” points are due to a crystal momentum 
fractionalization. We suggest neutron scattering mea¬ 
surements to detect low-energy intensity peaks around 
M" points to experimentally discover the phenomenon 
of crystal momentum fractionalization. This will be a 
smoking gun for the DSL in Herbertsmithite. 

We start with KAFHM for spins in Herbertsmithite 

= ( 1 ) 

(b> 

where J ~ 17 meV and the summation runs over nearest 
neighbor bonds. We will use the Schwinger fermion rep¬ 
resentation for spin-1/2 operator, Sf = 

Here are Pauli matrices. The fermionic spinon 

operator describes a spin physical Hilbert space 

within one-particle-per-site constraint flafia = 1- 

A spin liquid state is characterized by a mean-field 
Hamiltonian 

= ( 2 ) 

(iJ) 

The GPWFs for a spin liquid ground state and spin-1 
excited states are written as 

I'l')= 'Pci'i'U), = 'Pcfye.ii'i'U), (3) 

where Vq is the Gutzwiller projection operator to en¬ 
force one-particle-per-site constraint and |d>^p) is the 
mean field ground state, fe-a- is the operator for the 
wave packet with mean field energy level in the mean 
field Hamiltonian. 

Different choices of Xij in E*l- (2) give us different spin 
liquid states. FSL is a zero-flux state and has a large 
spinon Fermi surface. DSL has flux tt in the hexagons 
of kagome lattice. RSL has a random quenched gauge 
field aij on the bond, Xij = with —tt < Oy < tt 

randomly. 

From the ground-state GPWF, we can calculate the 
equal-time spin structure factor S'(q) 

= ( 4 ) 

ij 


where — Vj and the position summation runs over 

all sites on the kagome lattice. (■ ■ ■ )o is average over spin 
configurations in the ground-state GPWF. For a 12 x 
6x3 lattice, we specify the general periodic boundary 
conditions on the lattice 

h+L^ = fi, fi+Ly = ko = v^/2, (5) 

to get a full shell of mean field energy levels. The kagome 
lattice has the primitive basis ai 2 = -I- The 

reciprocal primitive vectors are gi _2 = ±27rka; -|- 
indicated by the purple parallelogram in Fig. 1 (d). S'(q) 
is periodic in extended MBZ (solid hexagons in Fig.l). 

In Fig. 1 (a), (b) and (c), we compare S'(q) among 
FSL, DSL and RSL. S'(q) has a similar overall feature 
for three spin liquid states. The main differences are the 
peak positions: While FSL and RSL have peaks around 
K points, DSL has peaks at M points in MBZ. As shown 
in Fig. 1 (e), along high symmetry directions in MBZ, 
FSL has a dip around M points and the dip feature is 



FIG. 1. Contour plots of equal-time spin structure factor S(q) 
as a function of momentum in FSL (a), DSL (b) and RSL (c). 
The solid and dashed hexagons are the magnetic extended 
and elementary Brillouin zones, respectively as shown in (d). 
Two spinon dispersion is periodic in the shadow parallelogram 
in (d). (e) S'(q) along high symmetry directions for FSL, DSL 
and RSL. 
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smeared out in RSL. DSL has a kink around M” points. 
RSL has no translational symmetry and S{q) is obtained 
as q-Fourier transformation in Eq. (4) on a 12 x 6 x 3 
lattice. 

Neutron scattering experiments measure the dynamic 
structure factor 

^ (T.5-(r)5+(0))o. 

° ij 

The projected Hamiltonian system within a subspace 
consisting of spin-1 excited states[46, 49] is given as 

= 0{i'f,ij) = {i'f\ij), (6) 

where \ij) is in Eq. (3). The matrix ele¬ 

ments in Eq. (6) are evaluated by using Monte Carlo 
methods[46, 49]. The projected Hamiltonian system 
{H, O} is diagonalized through a generalized eigen equa¬ 
tion, 

= enO\4>n) , (7) 

where \4'n) and e„ are spin-1 two-spinon wave functions 
and energy levels, respectively. In terms of them, we has 
the spectral representation 

5(q,a;) = ^ <5(a; - (e„ - (8) 
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FIG. 3. (a) Plot of integrated 5(q, oj) up to 0.6J and (b) 

plot of S(q, ca) with fixed frequency ujq = 0.03J, 0.18J, O.OIJ 
for FSL, DSL and RSL, as a function of momentum along 
high symmetry directions. 
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FIG. 2. (a) and (b) are contour plots of ^(qja;) in FSL (a) 
and DSL (b) as a function of frequency and momentum along 
high symmetry directions. The white solid circles are the 
lower edge i?edge(q) of the SWG. (c) and (d) are contour plots 
of S'(q, ta) for RSL with fixed energies ca = 0.05J (c) and 
ta = 0.5J (d) as a function of momentum. 


where Cq is the ground state variational energy[49]. 

FSL and DSL have translational symmetry and the 
projected Hamiltonian system {H, O} will be labeled ac¬ 
cording to the momentum in a relatively large system 
(12 X 6 X 3). The mean field dispersion for spinons has a 
finite-size gap on the 12 x 6 x 3 lattice with the bound¬ 
ary conditions in Eq. (5). The finite-size spin gaps are 
= 0.0291x1 and = 0.586lx| for FSL and DSL, 
respectively. Here jxj is the mean field spinon hopping 
amplitude. Correspondingly, spin-1 excitations have a 
gap, Eg = 0.03J and Eg = 0.18J for FSL and DSL, re¬ 
spectively. Contour plots of «S'(q, w) are shown in Fig. 
2 (a) FSL and (b) DSL, with broadening rj = 0.15J for 
delta function in Eq. (8), S{uj — e) —)■ The 

RSL state has no translational symmetry and the com¬ 
putation complexity increases considerably, ^(q, w) for 
RSL is computed in the whole Brillouin zone only on a 
4x4x3 lattice. In Fig.2, we plot S'(q, w) for RSL with 
fixed frequencies (c) uj = 0.05J and (d) w = 0.5J. 

The projected Hamiltonian system {H, O} for spin-1 
excited states has the SWC width around IFswc — 3J 
for three spin liquid states. In Fig. 3 (a), we plot the 
integrated S'(q, w) over 0 to 0.6J along high symmetry 
directions. Note that Fig. 1 (e) and Fig. 3 (a) has the 
same unit. The integrated S'(q, w) over 0 to 0.6J has 
about 20% intensity of the fully integrated of spin spec¬ 
tral weight <S'(q). To explore the low energy features, 
we plot 5'(q, w) with fixed energies a little higher above 
the finite-size gap, ujq = 0.03J, 0.18J, O.OIJ for FSL, DSL 
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FIG. 4. DSL: the lower edge of SWC, Sedge (q) — Eg, as a 
function of momentum along high symmetry directions. Eg 
is the finite-size gap. The dashed line is the lower edge of 
SWC from the meanfield theory with \x\ = 0.43J for hoping 
amplitude. 

and RSL, along high symmetry directions. The low en¬ 
ergy cones at M and M” points are very unique for DSL 
and clearly resolved in Fig. 3 (b). 

We can decompose S (q, w) into a spin matrix element 
M(q, w) and a density of two spinon excited states 

S'(q,w) = M(q,a;)i:)(q,w), (9) 

with ~ (en(q) - eo)), where e„(q) is 

the generalized eigen value of the projected Hamiltonian 
system {H, O} for a given momentum q. The lowest 
generalized eigen values ei (q) gives the lower edge of the 
SWC, Fledge(q) = ei(q) — Co, which is plotted as white 
solid circles in Fig. 2 (a) FSL and (b) DSL. 

FSL has a large spinon Fermi surface and its low energy 
sectors of spin-1 excited states strongly depend on the fi¬ 
nite size. In the thermodynamic limit K and K' points 
are equivalent under 60° rotation symmetry for the pro¬ 
jected Hamiltonian system {H, O}. However, in Fig. 2 
(a), the lower edge at K has a higher energy than K', 
Eedg,e{K) > Fledgesince the lattice shape (12 x 6 x 3) 
and the boundary conditions in Eq. (5) break the rota¬ 
tion symmetry. The lower edge for DSL resembles that 
in the mean field calculations. In Fig. 4, The lower edge 
of the SWC Fledge(q) — Eg for DSL fits well the mean 
field calculation with a finite-size spin gap Eg = 0.18J 
and the mean field hopping amplitude |x| = 0.43J. 

Different from FSL, fermionic spinons in DSL carry a 
crystal momentum fractionalization[6, 34]. Due to tt flux 
in the hexagons of kagome lattice, translational operators 
for spinons along the primitive lattice vectors ai .2 anti¬ 
commute with each other 

T 1 T 2 =—T 2 T 1 , Ti^2(x) = X-f ai.2. (10) 

As a result, the projected Hamiltonian system {H, O} 
has the spin-1 spectrum e„ with a period of one-quarter 
of the Brillouin zone. In other words, the momentum 
resolved density of states for spin-1 excitation continuum 
(ie. two-spinon excitations) has a period of one-quarter 


of the Brillouin zone. [4, 6]. Note that M' and M” points 
are equivalent to F points and cut the elementary BZ 
(shadow parallelogram in Fig. 1 (d)) into four pieces. 

The spin matrix element M(q, w) in Eq.(9) is periodic 
in the MBZ and M, M', M” and F are not equivalent any 
more in S'(q, w). While the magnetic intensity at F and 
M' points are suppressed in ^(q, w) as shown in Fig. 2, 
M and M” are still visible. The low energy intensity 
at M and M” in ^(q, w) is the implication of a crystal 
momentum fractionalization in DSL which is detectable 
in the neutron scattering measurements. 

Here we make several remarks on comparison between 
experiments and theoretical results. The three different 
spin liquids state have a similar overall shape in S'(q, w) 
with general features: a SWC spectrum over large en¬ 
ergy ^ 3J with low intensity in the elementary BZ and 
high intensity in 2nd BZ, in good agreement with exper¬ 
imental observations. In one-dimensional antiferromag¬ 
netic spin-1/2 chain, D[q,uj) is finite at the lower bound¬ 
ary and S{q,uj) has a divergent sharp lower edge due to 
M[q,uj)[AS]. Although enhanced at low energies, S'(q, w) 
does not diverge at the lower edge of the SWC in spin 
liquid states on the kagome lattice. So nearly invisible 
lower edge in Herbertsmithite experiment may not be a 
big issue. For FSL, S'(q, w) has a gap at the M point 
in the extended MBZ, in contrast to experiments. RSL 
also has low energy intensity peaks at K points inconsis¬ 
tent with experiments. In contrast, DSL has no obvious 
conflict with experimental observations. Due to a mo¬ 
mentum fractionalization, S'(q, w) of DSL has two low 
energy Dirac cones at M and M” points in the MBZ. In 
the experiments, below 1.5 meV, high intensity M points 
in ^(q, w) are connected through M” points instead of K 
point. So low energy intensity peaks at M” points due to 
a crystal momentum fractionalization may already be ob¬ 
served in experiments; however, these features are inter¬ 
preted as the impurity effects[19]. In the presence of im¬ 
purities, the system in Herbertsmithite may have low en¬ 
ergy gauge field fluctuations. We find that high intensity 
peaks at M and M" are stable against quenched gauge 
field fluctuations although the low energy boundary edge 
below SWC is smeared out [49]. Recently, Barlowite with 
AF ordering temperature Tn = IbK is studied as an¬ 
other kagome antiferromagnet[50]. Its non-magnetic (Mg 
or Zn) doped variety is proposed to has less imperfec¬ 
tions than Herbertsmithite[50, 51]. The new material is 
promising to clear the impurity issues. 

In conclusion, we study the fractional spin-wave con¬ 
tinuum in spin liquid states on the kagome lattice. We 
find out that DSL describes the experiments in Her¬ 
bertsmithite well. Besides a fractionalized spin moment, 
fermionic spinons in DSL carry a fractionalized crystal 
momentum which is also potentially detectable in fur¬ 
ther experiments. 
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I. SMEARED DSL 


In the Herbertsmithite compound, the system may not reach its true ground state due to the presence of impurities 
and the gauge fluctuations is quenched. We would like to study smeared DSL (SDSL) in this section. 

Xij has low-energy gauge fluctuations, Xij = i where a^- behaves as gauge bosons. The gauge fluctuations are 

very soft when the ground state of kagome spin model is close to a quantum phase transition. In the Herbertsmithite 
compound, the system may not reach its true ground state due to the presence of impurities and the gauge fluctuations 
aij is quenched. SDSL is described as 

l'I'SDSL)=^Gl'f'Sr^'), (1) 


where Xij is specified as the DSL and the gauge field —0.27r < < 0.27r is chosen randomly. For the SDSL, GPWFs 

for excited states are given as 




( 2 ) 


The 5'(q, w) is shown in Fig. 1 where the low energy intensity peaks at M and M” points are clear resovled. 


II. DETAILS FOR MC 


A. GPWFs for Spin liquids 


The mean held Hamiltonian in Eq. (??) can be easily diagonalized on the hnite lattice and it has N energy levels 
for both spin-up and down spinons, < •• • < Filling the lowest N/2 spin-up and N/2 spin-down 

energy levels, we obtain GPWF 


|G) = Vg (e 




t 

o ' 


2 _Jvac) G • • • ej^^^_^\vac] 


(3) 


which is assume as the ground state of the spin model (??) on the kagome lattice. This is the basic gound-state 
assumption. 

The ground GPWF |G) is a singlet state with Stot = 0. Since we are interested in the spin dynamics, we need 
construct excited states with total z-component spin = 1. Based on the mean held energy levels, we can construct 
many different states with total = 1 


/ 


= 1)^Vg 


t 


• eL/, Ivac) (g) e\ 


e] I vac) 

30 3N/2-2 ' ' 


(4) 


\ JV/2+1 7V/2-1 

where € {el, • • • , e]v_ il and 4 G {4g" ,4-i}- The dimension of sub Hilbert space with 5^^^. = 1 is 

2 


Dim(Hf,n[^(ot = l])~ ■ 


(5) 


^ r\ 
0 



FIG. 1. Smeared DSL: Contour plot of dynamic spin structure factor ^(q,^;) as a function of momentum at fixed frequency 
for (b)tu = 0.05J and (c) uj = 0.5J. 
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Excited-state assumption: the second assumption is that we are only interested in excited states which can be 
constructed by spinon particle-hole excited states 

l^tot = l)=^G(eI^e,;|vI/)) (6) 

where |^') = ej • ■ ■ ''' ejv/ 2 - 1 and G {e]v/ 2 ’ ’ ’ ’ >e]v_i}, G {ej^,-- ■ .e^/ 2 _i}- This truncated 

sub-Hilbert space has the dimension 

Dim(H[5t^t = l])~iVV4. (7) 

Note in Eqs. (5) and (7) we use ^ instead of =, because those states may be not linearly independent. The dimension 
in Eq. (5) ~ 0(4:^), N —)■ 00 is much larger than the full spin Hilbert space dimension 2^ when N is large. On 
square lattice for uniform RVB state, only around half of the excited states in Eq. (6) are linearly independent and 
Dim(H[5t^,t = 1]) ~ iVV8. 


B. Dynamic spin strncture factor 

At zero temperature, the spin-spin correlation is given as 

{TrS-iT)S+iO))o = {G\e^^S-e-^^S+\G). ( 8 ) 

Under the above two GPWE assumptions, we can expand the Hamiltonian in the truncated Hilbert space 

H = J2En\En){E„\, (9) 

n 

where \En) is the normalized orthogonal basis. Then we have the dynamical spin structure factor 

^(q,.;) = - (£!„ - £;o))|(E„|5+|G)|2, (10) 

n 

where the translational symmetry is assumed. 


C. Projected Hamiltonian system in truncated sub-Hilbert space 

We are working on the truncated sub-Hilbert state with total spin = 1 in Eq. (6). We should project the spin 
Hamiltonian onto the truncated Hilbert space to obtain the expansion as in Eq. (9). 

Due to projection, the Gutzwiller-projected particle-hole states in Eq. (6) are not orthogonal any more (even not 
linearly independent). To obtain the normal orthogonal basis, we need calculate the overlaps between these states 
and the Hamiltonian matrix 


0 (bj) = (*|j). = {i\H\j), (11) 

Diagonalize the overlap matrix 

OU(”)=A„U(”\ (12) 

we obtain the normal orthogonal basis 

N^/4-1 

Wn) = 

i=0 

Note here zero eigenvalues A„ imply the |a„) is not an independent state and should be removed from the orthogonal 
basis. The rank r of overlap matrix O is the dimension of the truncated sub-Hilbert space, Dim('H[S'fot = 1]) = r. 
Generally, r <= iV^/4, e.g., uniform RVB state on the square lattice has r ~ N'^/8. 

Based on the orthogonal basis \an), we can obtain the Hamiltonian 


a/ An 


^)V, 


(n) 


(13) 




( 14 ) 
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Diagonalize the above Hamiltonian matrix, we can obtain the normal orthogonal expansion of Hamiltonian in Eq. (9). 
When the overlap matrix O has the full rank, the above procedure can be re-expressed in the generalized eigenvalue 
problem^ 

= £^^4/^ ( 15 ) 

where |i?„) = Therefore the dynamic spin susceptibility is given as 

- (£;„ - Eo))| P, 

n ij 

( 16 ) 

where the spin operator is expressed as 5+ = Wv’f • 

For reasonably large system size, e.g. = 12 x 6 x 3 on the kagome lattice, the dimension Dim('H[S'tot = 1]) is 
still huge for numerical calculations and large memories are needed to store O and El matrices. If the system has the 
translational symmetry, we can decompose Dim('H[S'fQj, = 1]) according to the lattice momentum q 

= 1 ]) = 0 Bim{n[S4,, = l; q]). ( 17 ) 

q 

Particularly, for FSL state, Dim(77[S'fot = 1; q]) are not the same for different momentum q. The typical value of the 
dimension is Dim(77[S'fot = 1; l]) ^ 1^(-^)) which is small enough for the efficient numerical calculations. 


D. Monte Carlo algorithm 


The key evaluations are the matrix elements for O and El matrices in Eq. (11). This can be done by Monte Carlo 
method. We will take the sampling strategy developed by Li and Yang in Ref. ? : using single Markov chain to 
generate spin configurations for all element evaluations: 


o(bj) = 

a 


p{a) 

{i\a){Ha\j) 


p{a) 


P(a). 


(18) 


where |a) is spatial spin conhgurations generated according to the Monte Carlo sampling probability p(pi). Given the 
spatial spin configuration |a), Gutzwiller-projected particle-hole states in Eq. (6) are all Slater determinants 

(a|i) = (a^|t^) X (a^lt^). (19) 

For the evaluations, we would pick up one reference state which is the lowest mean field particle-hole state 

\R) = "Pg (ej • • • e]v/2|vac) O . (20) 

and calculate its Slater determinant 

{a\R) = (at|i?t) X (21) 

Since every determinant differs with only by one column, it is easy to evaluate using rank-1 

determinant update. The computation complexity is 0{N). The Hamiltonian matrix element determinant is 

{Hoi^/\\i\/\) (22) 

which differs with by one row and one column and also can be evaluated using rank-2 determinant update. 

The computation complexity is 0{N^). So the full problem has the complexity is 0{N^). 

In Refs. ? , the total weight of all states with the same momentum q is used for /9q(a) = l(“l*)q)P- During 

the update according the total weight, the matrices for the reference state determinants (a^|i?^) or (a^|i?^) may be 
singular. To avoid the singularity, we use the probability function 

p{a) = Katl-Rt) X (a^|i?^)|A 


(23) 



